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Using the first law of binary black-hole mechanics, we compute the binding energy E and total angular 
momentum J of two non-spinning compact objects moving on circular orbits with frequency £2, at leading 
order beyond the test-particle approximation. By minimizing E(£i) we recover the exact frequency shift of the 
Schwarzschild innermost stable circular orbit induced by the conservative piece of the gravitational self-force. 
Comparing our results for the coordinate invariant relation E(J) to those recently obtained from numerical 
simulations of comparable-mass non-spinning black-hole binaries, we find a remarkably good agreement, even 
in the strong-field regime. Our findings confirm that the domain of validity of perturbative calculations may 
extend well beyond the extreme mass-ratio limit. 
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Introduction. — The problem of motion has always played 
a central role in physics, especially that of binary systems in 
gravitational physics. While in Newtonian gravity the two- 
body orbital motion can be solved analytically, the exact gen- 
eral relativistic solution for two black holes was only obtained 
very recently, by means of numerical simulations (T). 

The post-Newtonian (PN) approximation to General Rel- 
ativity has proven a valuable tool to describe the dynamics 
of widely separated inspiralling compact-object binaries El , 
which are expected to be the main gravitational-wave (GW) 
sources for existing ground-based interferometric detectors 
(LIGO/Virgo) and future space-based antennas. Furthermore, 
it has achieved successes ranging from Solar-system tests of 
gravity [3], to measurements of the damping rate of binary 
pulsars due to GW emission H, and to remarkable practical 
applications such as the Global Positioning System (GPS) [5|. 
However, the PN approximation becomes inaccurate during 
the late stages of the binary's inspiral, and breaks down dur- 
ing the final plunge and merger. While numerical-relativity 
(NR) simulations can describe these highly relativistic phases 
of the evolution, they are still too time-consuming to construct 
GW template banks covering the whole parameter space of 
compact binaries. 

An alternative analytical approach that can improve our 
knowledge of this highly relativistic regime is the gravitational 
self-force (GSF) formalism, a natural extension of black-hole 
perturbation theory |6|. The GSF approach relies on an ex- 
pansion in the binary's mass ratio, and is the natural tool to 
model extreme mass-ratio compact binaries, which are among 
the most promising GW sources for space-based detectors J7). 
Unlike the PN approximation, which breaks down when the 
binary's velocity becomes close to the speed of light, the GSF 
formalism remains valid for highly relativistic systems, even 
in the strong-field regime. Moreover, it has recently become 
clear that the GSF can also describe at least certain aspects 
of the dynamics of comparable-mass binary systems, such as 
the relativistic periastron advance [8|, possibly allowing the 
construction of vaster and more accurate template banks that 
would be crucial for GW astronomy. 



In this Letter, we confirm that picture by computing the 
binding energy E and angular momentum J of a non-spinning 
circular-orbit compact binary system within the GSF formal- 
ism. Our results provide a surprisingly accurate description 
of comparable-mass binaries, as we verified by comparing the 
coordinate-invariant relation E(J) that we obtain with recent 
NR data J9). In addition, they pave the way to adiabatic evo- 
lutions of extreme mass-ratio inspirals that include the effect 
of the first-order conservative GSF. 

More specifically, we start from the first law of mechanics 
for binaries of spinless compact objects moving along circu- 
lar orbits, and modelled as point particles. This relation, re- 
cently established in Ref. [ 10], gives the variations of the total 
Arnowitt-Deser-Misner (ADM) mass M and total angular mo- 
mentum J of the binary system in response to small variations 
of the individual masses ntA (A = 1,2) of the compact objects 
according to 

8M-Q.8J = zi8m l +Z2dm 2 . (1) 

Here we set G = c = 1 (a choice that we adopt throughout 
this Letter), £2 is the circular-orbit frequency, and za are the 
so-called "redshift observables", namely the gravitational red- 
shifts of light rays emitted from the particles, and received far 
away from the binary along the direction perpendicular to the 
orbital plane ifTTI . In a convenient gauge, the redshifts sim- 
ply coincide with the inverse time components of the four- 
velocities u" of the particles, namely za = 1/ m a ATI . The 
relation ([T} is the point-particle analog of the celebrated first 
law of black-hole mechanics 8M — £2h 8J = 4K"m UT 5to ut lfl2l . 
where m m = yjA/{\6K) is the irreducible (or Christodoulou) 
mass of a black hole of mass M, spin J = aM, surface area A, 
uniform surface gravity fc, and horizon frequency £2h. (See 
Refs. iflOl [T3l for more details on the first law of black-hole 
mechanics, or "thermodynamics".) 

In this Letter we show that the first law ([1} can be used, 
in conjunction with existing perturbative GSF calculations of 
the redshift observable, to compute the binding energy E = 
M — (m\ + m%) and angular momentum J of compact binary 
systems on circular orbits, at leading-order beyond the test- 
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particle approximation. (See Ref. [10] for a discussion of the 
applicability of Eq. ([1} to GSF calculations.) 

As an immediate application of our results, we recover 
the exact frequency shift of the Schwarzschild innermost sta- 
ble circular orbit (ISCO) induced by the conservative piece 
of the GSF, as computed previously from a stability analy- 
sis of slightly eccentric orbits around a non-spinning black 
hole lfl4l . We then compare our newly derived GSF-accurate 
expression for the invariant relation E(J) with the results 
of spinless binary black-hole simulations with mass ratios 
q = m\/m.2 = 1, 1/2, and 1 /3 (we assume m\ ^ mi through- 
out this Letter). Finally, we summarize our results and discuss 
future work. 

Binding energy and angular momentum. — The ADM mass 
M, total angular momentum 7, and redshift observables za are 
all functions of the circular-orbit frequency Q. and the indi- 
vidual masses iua- The first law ([T| thus implies dM/dQ. = 
QdJ/dQ and dM/dmA — QdJ/dniA = Za- Applying the 
change of variables (Q.,mi,nt2) — > (m,V,x), where m = mi + 
W2 is the total mass, v = mimi/m 2 = pt/m the symmetric 
mass ratio, and x = (mil) 2 / 3 the usual dimensionless invariant 
PN parameter, these relations can be combined to give 
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where j% = M — £27 can heuristically be viewed as the energy 
of the binary in a co-rotating frame, and A = {mj — nt\)/m = 

fEJ, 



y/T— 4v is the reduced mass difference. To derive Eq. 
we used the fact that the dimensionless ratio Jm is a func- 
tion of v and x only (cf. the discussion following Eq. (4.21) of 
Ref. ifTUl ). A similar equation for particle 2 can be obtained 
by substituting A — > —A. 

At first order beyond the test-particle approximation, the 
following mass-ratio expansions hold for the redshift zi, the 
specific binding energy E = (M — m) / jj., the dimensionless 
angular momentum 7 = 7/ (mjl), and M = (./# —m)j\i: 



z\ = V1-3x + vzsfW + ^(v 2 ), 
( \-2x 



V\/l -3x 
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y/x(l-3x) 
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The lowest-order terms are the well-known Schwarzschild re 
suits in the test-mass limit. Since v = q/(l + q) 2 = q -\ 
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, these equations remain valid if the symmetric mass 
ratio V is replaced by the usual mass ratio q. Substitut- 
ing the expansions |3j in Eqs. 0, one obtains the follow- 
ing relations between the various GSF corrections: zsf(*) = 
2^f SF (jc) - 2(y1 -3x- l)-x/y/l-3x, E SF (x) = Jt<x(x) - 
T-^sf( x )' anc l ^sf(^) = — jtj-^sfM' wnere we denote 
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The knowledge of the GSF correction zsf to the redshift (and 
thus of its first derivative Zg F = dzsp/dx) therefore immedi- 
ately gives the corrections £sf an d -^sf to the test-particle re- 
sults for the binding energy and angular momentum. 

The GSF effect on the coordinate-invariant relation zi(x) 
was first computed in the Regge- Wheeler gauge in Ref. ifTTl . 
Alternative GSF calculations based on different gauges (har- 
monic gauge and radiative gauge) were later found in agree- 
ment to within the numerical uncertainties |fT31 [T6ll . Collect- 
ing all the published GSF data IfTTl fTsHrTl . we have at our 
disposal 55 data points for zsf(x), with x ranging from to 
1 /5, and with relative errors lower than 10~ 6 . In particular, 14 
of these points lie in the "strong-field interval" 5m ^ r ^ 10m, 
where r = m/x is an invariant measure of the separation. 

Following Ref. Ifl8l . it is convenient to represent these 
data by a compact analytic expression, in terms of a ratio 
of polynomials in the PN parameter x. We adopt the model 
Zsf(x) = 2x(l + a\x + a2X 2 ) / {\ +aix + a4X 2 +asx 3 ), which 
accounts for the asymptotic behavior zsf(x) = 2x + (?(x 2 ) 
when x — > IfTUl . Performing a standard least-squares fit, 
we find that the coefficients a\ = —2.18522, ai = 1.05185, 
a 3 = -2.43395, a 4 = 0.400665, and a 5 = -5.9991 reproduce 
the data to within 10~ 5 . We notice that the extrapolation of our 
fit beyond the data point with the smallest separation (r = 5m) 
diverges very close to r = 3m. This suggests that the GSF con- 
tribution to z\ may have a pole near the Schwarzschild circular 
photon-orbit (or "light-ring"), located at xlr = 1/3. 

Using this fit for zsf(*) in Eqs. (HJ, the GSF contributions 
Esf and 7sf to the binding energy E and angular momentum 
7 can easily be computed in the range ^ x ^ 1/5, with a 
comparable accuracy. 

Self-force correction to the Schwarzschild ISCO. — We re- 
call that the Schwarzschild ISCO is defined by the onset of a 
dynamical (radial) instability for circular orbits. As a first ap- 
plication of our results, we show how the GSF-induced shift 
of the ISCO frequency can be recovered very simply from the 
expressions pb) and ( f4a[ > for E(£i). 

Beyond the test-particle approximation, the orbital fre- 
quency of the ISCO is given by 



mOisco = 6~ 3/2 [1 + vC n + ^(v 2 )] 



(5) 



where the coefficient Cn, which encodes the effect of the 
conservative piece of the GSF, has recently been computed 
by Barack and Sago (BS) (HI [19). Performing a stability 
analysis of slightly eccentric orbits near r — 6m, they found 
C^ s = 1.2512(4) (19). This strong-field benchmark has since 
then been used as a reference point for comparison with other 
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analytical and numerical methods 11201 1211 . and for calibra- 
tions of the effective-one-body (EOB) model lfl8l I22ti24l . 

On the other hand, the minimum-energy circular orbit 
(MECO) is defined as the minimum of the binding energy 
E(£l); the MECO's orbital frequency iijyiECO thus satisfies 



dE 
da 



^MECO 



= 0. 



(6) 



It was shown in Ref. [25 1 that the notions of ISCO and MECO 
are formally equivalent; hence £2isco = ^meco- This result 
does not rely on any PN expansion or perturbative analysis, 
and thus holds for any mass ratio, even in the strong-field 
regime; it only requires that the binary's dynamics can be de- 
rived from a Hamiltonian. Inserting Eqs. ( |3b| > and ( |4a| i for the 
binding energy E(£l) in the condition d§), we then find the 
following expression for the ISCO frequency shift: 
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In order to compute the GSF correction to the Schwarzschild 
ISCO, one thus only needs the first and second derivatives of 
the GSF correction to the redshift, evaluated at x — 1 /6. By fit- 
ting our full data set for zsf (x) (or the strong-field subset 5m ^ 
r ^ 10m only) to different models, we find Ca = 1.2510(2), 
which agrees with BS's result C^ s at the Iff level. 

Moreover, Eq. |7]) allows for a highly accurate determina- 
tion of the ISCO frequency shift. Indeed, while previous cal- 
culations relied on a stability analysis near a singular point, 
using a GSF code capable of handling slightly eccentric or- 
bits, Eq. |7]) only requires the evaluation of a regular function 
near x = 1/6, using a much simpler GSF code for circular 
orbits. Current circular-orbit GSF codes implemented in the 
frequency domain can already deliver highly accurate results: 
for instance, the data reported in Refs. ifTTl l26l for the GSF 
correction to u\ = l/z\ (for separations r ^ 200m) are accu- 
rate to within 10~ 13 , and similar accuracies should be achiev- 
able at least down to r = 5m. With high-accuracy GSF data 
for zsf near x = 1/6, it will become possible to determine the 
ISCO frequency shift Cq much more accurately than ever be- 
fore; a valuable result given the physical significance of this 
genuinely strong-field effect. In particular, this could prove 
useful to cross-check the results of different GSF codes. 

Comparison with numerical relativity. — We now combine 
our results ([3]l and Q for the binding energy E(£l) and angu- 
lar momentum J(£l) to compute the coordinate-invariant rela- 
tion E(J), at leading-order beyond the test-mass approxima- 
tion, and compare it to the results recently obtained in Ref. [9 1 
using accurate NR simulations [27 1 (with Cauchy characteris- 
tic extraction [28 1) of non-spinning black-hole binaries with 
mass ratios q = 1,1/2, and 1/3. This comparison is pre- 
sented in Fig. [TJfor the q = 1 case, and gives similar results 
for q = 1 /2 and q=l/3. 

The curve labelled "Schw" shows the relation E(J) for a 
test mass on an adiabatic sequence of circular orbits around 
a Schwarzschild black hole, given in parametric form by 
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FIG. 1: In the upper panel, the (specific) binding energy E = Ejp. of an 
equal-mass black-hole binary is shown as a function of the (dimensionless) 
total angular momentum / = as computed in numerical relativity 

("NR"), in PN theory ("3PN"), in the EOB model ["EOB(3PN)"], in the test- 
particle approximation ("Schw"), and including the conservative gravitational 
self-force ("GSF<f and "GSFv"). The 3PN, EOB(3PN), and test-particle 
curves show cusps at their respective ISCO/MECO; the lower branches corre- 
spond to stable circular orbits, while the upper branches correspond to unsta- 
ble circular orbits. We might find similar branches for the GSF curves when 
more data closer to the light-ring become available. The "GSFg" and "GSFv" 
curves are only shown in the region where numerical data for the self-force 
are available (i.e. r ^ 5m, corresponding to /~ 3.126 for the "GSFv" model 
and to J ~ 1.896 for the "GSFg" model). The differences between the vari- 
ous models and the NR result, AE = E — £nr, are shown in the lower panel 
(down to the ISCO/MECO, when that is present), with the exception of the 
"GSFg" model, which quickly grows beyond the plot range. The shaded area 
represents the error affecting the NR results. 



the lowest-order terms in Eqs. ( |3b] i and ( |3c| l. It necessar- 
ily presents a cusp at the Schwarzschild ISCO, located at 
m 

Qfchw _ 068 As 

can be seen in Fig. [I] this result is in 
reasonable agreement with the NR data al the way down to 
the ISCO, where the difference reaches AE ~-1.7x 10~ 3 . 
The astonishing agreement between the test-mass result and 
the NR result for q = 1 suggests that the relation E(J) should 
be almost independent of q, at least for "large" orbital sepa- 
rations. Indeed, we verified this using the NR data for q = 1, 
1 /2, and 1 /3, as well as the other models shown in Fig. |T| 

The E(J) adiabatic relation given in parametric form by 
Eqs. ( |3b| i and ( |3"c] > includes the effect of the conservative GSF. 
Expressing Eqs. ([3b} and ( |3c] i in terms of the mass ratio q 
(which simply amounts to replacing v — > q), we obtain the 
curve "GSFg". Given the large mass ratio involved (q = 1), 
the poor agreement with NR is expected, and the agreement 
does not improve significantly forg = l/2org = 1/3. How- 
ever, the GSF result expressed in terms of the symmetric mass 
ratio v ("GSFv") compares remarkably well with the NR re- 
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suit, with a difference that grows larger than the numerical 
error only near r = 5m (which corresponds to J ~ 3.126 for 
the "GSFv" model, where it reaches AE ~ 3.5 x 10~ 4 ). 

For completeness, we also show the invariant relation E(J) 
as computed in the adiabatic PN approximation ("3PN") and 
in the EOB adiabatic model ["EOB(3PN)"]. The PN re- 
sult is given in the parametric form {.c(fl), /(£!)} by, e.g., 
Eqs. (5.13) and (5.8) of Ref. Q9\, where we set V = 1 /4, and 
use the known values ©static = and ©kinetic = 41 /24 for the 
3PN static and kinetic "ambiguity parameters". We observe 
very good agreement with the NR data all the way down to 
the cusp occurring at the 3PN MECO (mCl™ co ~ 0.129), 
where the difference grows to AE ~ — 10~ 3 . In that respect, 
we point out that this PN approximant performs much better 
than the PN approximant considered in Ref. [9| [see Eq. (5) 
there], which is obtained through the series reversion E[Cl(J)] 
and subsequent re-expansion through 3PN order. Although 
technically correct, the resulting (uni-valued) function nec- 
essarily fails to capture the cusp at the MECO. Moreover, 
we checked that even in the test-particle limit the paramet- 
ric form {E(Cl),f(Cl)} is closer than the expression E[Cl(J)] 
to the Schwarzschild result. The poor behavior of some PN 
approximants, such as E[Cl(J)], is not a surprise in PN theory 
(see, e.g., Ref. Il30l ). The EOB result in Fig. |T] is produced 
using the 3PN model of Ref. [31], with a (1,3) Pade model 
for the effective metric component gff (see, e.g., Ref. J9]). As 
can be seen the difference with respect to the NR data grows 
as large as AE ~ —9 x 10~ 4 near the cusp occurring at the 
EOB ISCO (mflgg, ~ 0.088). 

We emphasize that the NR curve was obtained by Ref. ||9) 
from an actual binary black-hole evolution, and therefore 
includes non-adiabatic effects during the late inspiral and 
plunge. These effects are not captured by our adiabatic mod- 
els, and may in part explain the differences from the NR result 
in Fig. [I]at small J (i.e. at large CI). We also note that the re- 
markable agreement between the adiabatic models and the NR 
data for E(J) does not automatically imply that the same will 
hold true for the invariant relations E(Cl) and J (CI). 

Finally, we stress that although the GSF-accurate binding 
energy necessarily has an ISCO/MECO for small mass ra- 
tios (since for q = it reduces to the binding energy of a test 
mass in a Schwarzschild background), it does not present an 
ISCO/MECO (at least for r > 5m) for q = 1, 1 /2, and 1/3. It 
remains to be seen if this holds true even when GSF data for 
r < 5m become available. 

Summary and future work. — Recently, the general relativis- 
tic periastron advance of non-spinning black-hole binaries on 
quasi-circular orbits was computed in NR, and compared to 
the prediction of the GSF (as well as to other approximation 
techniques) [8]. By expressing the GSF result in terms of the 
symmetric mass ratio v rather than the usual mass ratio q, the 
GSF prediction was found in remarkable agreement with the 
exact NR result, even for comparable masses. This prompted 
the authors of Ref. [8 1 to suggest that the domain of validity of 
perturbative calculations may extend well beyond the extreme 
mass-ratio limit. Our new, alternative comparison based on 



the invariant relation E(J) strongly supports this expectation. 
A similar observation was previously made for the dissipative 
component of the GSF, based on a comparison of perturbative 
and NR calculations of the GW energy flux for head-on colli- 
sions 1 32]. The "scaling-up" procedure q — > v has also been 
used in the context of perturbative calculations of the linear 
momentum flux for quasi-circular orbits ll33l l34l . In the fu- 
ture, our analysis should be revisited using more GSF data in 
the very strong-field region 3m < r ^ 5m, as well as including 
dissipative effects in a consistent GSF evolution l35ll . 

Our expression for the binding energy E(Cl) can also be 
used to compute the EOB effective metric component gff ex- 
actly, through linear order in the mass ratio. Furthermore, by 
combining this result with the recent GSF/EOB comparison of 
Ref. 1T81 for the periastron advance in quasi-circular compact 
binaries, the g~f component of the EOB effective metric can 
also be computed. These results, which completely determine 
the EOB metric for spinless binaries through linear order in 
the mass ratio, are presented in the companion paper [36|. 

In the Schwarzschild spacetime, circular orbits for massive 
particles exist for any radius r > 3m. The redshift observable 
can thus be calculated, at least in principle, at any such radius. 
When more data for zsf(^) near x = 1/3 become available, 
our formulas ([3]) and Q for the binding energy and angu- 
lar momentum will provide information about the shift of the 
light-ring frequency induced by the conservative GSF acting 
on ultra-relativistic particles, or photons. 

In this respect, we emphasize that the connections estab- 
lished by Eqs. Q are particularly useful to explore the highly 
relativistic regime: while standard perturbative analyses can- 
not describe the binary's dynamics beyond the ISCO 11371 . 
the relations Q give direct access to the binding energy 
E and angular momentum / in the very strong-field regime 
3m < r ^ 6m, using only "routine" GSF calculations of the 
redshift z\ for circular orbits. 
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